Using clusters in SZE + X-Ray surveys 
as an ensemble of rulers to constrain cosmology 
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Ongoing and upcoming surveys in x-rays and SZE are expected to jointly detect many clusters 
due to the large overlap in sky coverage. We show that, these clusters can be used as an ensemble 
of rulers to estimate the angular diameter distance, dA{z) . This comes at no extra observational 
cost, as these clusters form a subset of a much larger sample, assembled to build cluster number 
counts ^ . On using this dA{z) , the dark energy constraints can be improved by factors of 1.5 - 
4, over those from just ^ . Even in the presence of a mass follow-up of 100 clusters (done for mass 
calibration), the dark energy constraints can be further tightened by factors of 2 - 3 . Adding dA{z) 
from clusters is similar to adding dhiz), from the SNe observations; for eg., ^ (from ACT/SPT) 
plus dAiz) is comparable to ^ plus di^{z) in constraining Q,m and ag. 



PACS numbers: 95.35.-l-x, 98.80.-k, 98.65.Cw, 98.62.Py 



I. INTRODUCTION 

Large cluster surveys like the SPT, ACT, Planck and 
eROSITA promise to detect from a thousand to a few 
hundred thousand clusters in the coming decade. The 
abundance and redshift distribution 4r^ of these clusters 

az 

are important probes to understand the nature of dark 
energy as well as to constrain other cosmological param- 
eters like n„, and cts S i, H, IH [ll] . To deduce a cos- 
mology from these ^ observations, one requires a precise 
knowledge of the limiting mass of the survey as a function 
of redshift. One frequently uses proxy observables such as 
X-Ray surface brightness and temperature [6|, Sunyaev 
Zel'dovich effect (SZE) decrement [3,[2J1, cluster richness 
[l3,[20[ and lensing [23,|31| for the masses of clusters, re- 
lated through simple power-law scaling relations. These 
scaling parameters are highly degenerate with the cos- 
mological parameters, and breaking this degeneracy is 
crucial to obtain tight constraints on cosmology. This 
may realized, for example, through the so called 'self- 
calibration' techniques [^, [H, [l3| • Other approaches in- 
clude an 'unbiased' mass follow-up of a sub-sample of the 
survey clusters [16|, |T3| or better theoretical modeling of 
clusters topredict the form of mass-observable scaling re- 
lation 0, [22, [sol ■ O'^^ can also try to optimize the cluster 
surveys so as to get the best possible survey yield 0, [HI . 

Measurement of the angular diameter distance, (Ia, at 
the redshift of the cluster using a combination of SZE and 
X-Ray observations have been routinely made over the 
last 30 years. The results have suffered, in the past, from 
various systematics and reliable estimates have only been 
achieved recently with analysis of statistically significant 
samples of galaxy clusters [J, [2l| . These new observations 
have demonstrated the power of using clusters to mea- 
sure 'c^A vs z' and use it to study the expansion history 
of our Universe [la|. However, these recent progress has 



been done with targeted observations. Since, targeted 
observations are costly, this approach limits the size of 
the sample of < 100. In this letter, we show that one can 
build up an 'ensemble' of dA^z) by picking a sub-sample 
of clusters discovered in both X-Ray and SZE surveys 
with overlapping area and redshift coverage. Since, the 
surveys are already geared towards getting clusters for 
^, we get the dA{z) without any extra targeted obser- 
vations. Addition of dA{z) to ^ helps in tightening 
cosmological constraints, especially on dark energy equa- 
tion of state. This is not surprising, since using dA{z) 
from clusters is akin to adding di,{z) information from 
supernovae (SNe) observations. 

The rest of this paper is organized as follows. In § 2, we 
briefly discuss cluster number counts and the procedure 
to estimate dA{z) from SZE + X-Ray observations; in § 3 
we describe the surveys, the choice of cosmological and 
cluster models and also summarize our methodology; in 
§ 4, we forecast constraints on cosmological parameters; 
and finally, we conclude in § 5. 



II. PRELIMINARIES 



The redshift distribution of detectable clusters is given 
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by, 

dN , ^ ^r. dV , ^ f°° dn(M,z) ^,^^ ,,,^ ,^^ 

where dV/dzdfl is the co-moving volume element, 
f{M, z) is the cluster selection function incorporating 
a logarithmic scatter in the luminosity to mass con- 
version and -^ is the halo mass function taken from 
simulations 10). For a X-Ray/SZE survey the lim- 
iting mass Miii„(2:) is found from flux limit, /nm, of 
the survey. For X-Ray survey, we adopt luminosity- 
mass relations from Vikhlinin et al. [26| given by 

Lx - Ax{^Y''E^-^^{z){J^y'-''\l + zV-. For 
an SZE survey, our SZE flux-mass relation is Yd\ ~ 
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FIG. 1: (a) Miim as a function of redshift plotted for various X-Ray and SZE surveys. The dA{z) mock catalog is created 
by using the effective Miim(z), which is the maximum of the Mum corresponding to a given Nhcum, the Mum for determining 
temperature of the clusters detected through eROSITA and the Mum for detecting clusters through Planck/SPT. (b) A mock 
dA{z) catalog of 578 clusters (as found from an overlap of ACT/SPT with eROSITA) with 25% Gaussian errors plus scatter. 
The thick line indicates the cJa from fiducial cosmology. 



^szE ( a.M^V ) E^'iz) (1 + zy^^ , where Y is the 
integrated SZE distortion and d^iz) is in Mpc. The 
parameters 'yx/lszE mimics any 'non-standard' evolu- 
tion of cluster scaling relations [1^ . For the X-Ray sur- 
vey, our fiducial parameter values are: log(Ax) = —4.25, 
ax = 1.61 and ^x — with a log-normal scatter of 0.246; 
while the corresponding values for the SZE scaling rela- 
tion are: log(AszE) = 1-75, aszE = 1-61, ^sze = with 
a log-normal scatter of 0.2. 



The distance dj^{z), of a cluster observed at a redshift 
z depends on the expansion history of the Universe as, 
dxiz) = TjT^ Jq h(z'') '^lisre U(^z) is the Hubble expan- 
sion. The seminal paper by [2j| defines the procedure to 
measure dK{z) from observations of clusters in SZE + X- 
Ray. This depends on 'assuming' a profile for the Intra- 
cluster medium, ICM, typically taken to be an isother- 
mal /3-profile 32] with n^{r) = n^^[\ -t- (r/rc)^)"^''/^, 
where n^ is the electron number density, r is the radius 
from the cluster's center, Tc is the core radius, and /3 
is a power-law index. Note, that one can always use a 
more realistic ICM model with better data. For example, 
one can fit a temperature profile to improve upon the 
errors. The X-Ray brightness, Sx oc J ri^Aeedl, where 
Aee is the X-Ray cooling function of the gas; while for 
SZE, AT ex JueTedl, where Te is the ICM tempera- 
ture. Eliminating Ue, gives the angular diameter distance 

dA oc g'^™g'"' ; see Birkinshaw Q for more details. 



III. ESTIMATING CONSTRAINTS FROM 
FUTURE SURVEYS 



Fiducial cosmology, priors and survey 
descriptions 



We adopt our fiducial cosmology from the WMAP 
7- year results (Table 6 of [IJ]) along with the follow- 
ing priors (Ans,Arjfc,A/i) = (0.015,0.0037,0.028). For 
simplicity, we choose a flat Universe since for an open 
wCDM model, WMAP7-FBAO-FH0 tightly constraints 
Al^tot < 0.007. With a small number of clusters, as sug- 
gested by recent SZ observations 25], just ^ data would 
not be sufficient to break the cosmology-cluster physics 
degeneracies. We therefore put priors on the scaling pa- 
rameters AA = 0.003 and Aa — 0.015 motivated by 
recent observations [1| . In addition, we put a weak prior 
of A7 = 0.2 in all the cases. 

We consider the following surveys with overlapping sky 
coverage - 

(1) ACT/SPT: We model the ongoing ACT/SPT survey 
as a 4000 deg^ with a /nm of 75 mJy (at 150 GHz), so as 
to give < 1000 clusters. 

(2) Planck: Ongoing all sky SZE survey. We take fum 
to be 300 mJy (at 353 GHz) which returns < 2000 clus- 
ters in ~ 32000 deg^. The higher flux limit means that 
Planck would detect only massive low z clusters. 

(3) eROSITA: Upcoming full sky X-Ray survey. We as- 
sume a [0.5-2.0 keV] band /nm = 4 x 10~^^ erg cm"^ s"^ 
which gives us < 1 x 10^ clusters for ^^ 32000 deg^. All 
clusters detected by both Planck as well as ACT/SPT are 
expected to be detected by eROSITA as it has a much 



(see Fig. 1(a)) 



B. Methodology 



The redshift distribution of clusters is obtained by us- 
ing the cluster scaling relations to convert flux to the cor- 
responding lowest observable cluster mass Miini(z). Fig- 
ure 1(a) shows the Mum for some of the surveys that 
we have considered. We further place a lower cut-off of 
1.3 X 10"/i~iMq(M5oo) on Mum and a redshift cut-off 
of 0.1 [H. The ^ (Az = 0.1) likelihoods are computed 
using Cash C statistics. 

We generate the mock c^a catalog as follows. For a clus- 
ter to be visible in both X-Ray and SZE surveys, its mass 
must lie above the highest of the Miim(z) from either of 
the surveys. To estimate cIa, one also needs at least a 
single isothermal temperature measurement of the ICM, 
which is possible for 10 times the detectable flux; for this 
we calculate the corresponding higher Mii,„. Next, to fit 
a /3-model, we also require that the clusters be well re- 
solved; i.e. be larger than a certain minimum angular 
size, so as to estimate its core radius dc- This is fulfilled 
by the condition that 9c be at least A'boam times the min- 
imum resolution of the survey (16" for eROSITA) |34| . 
We convert this to isophotal size Rj = OjdA and use the 
scaling relation between cluster mass and its size [13 to 
enforce this constraint. 



Fig. 



1(a) shows all the limiting 
The 



masses described above as a function of redshift 
redshift c^a catalog is constructed from [1] by integrat- 
ing over the highest of these Miim(2)'s in the plot. The 
dA^z) from these clusters are distributed randomly with 
a Gaussian scatter of 25% about the dA{z) from fidu- 
cial cosmology (see Fig. 1(b) I. Such catalogs are created 



for overlaps of SZE and X-Ray surveys like ACT/SPT 
+ eROSITA and Planck + eROSITA for TVbeam's of 2, 
3 and 4. A higher value of A^boam implies a selection of 
only the larger clusters; see Fig. 1(a) The ^a catalog is 



analyzed using a chi-square statistic for the likelihoods. 

Finally, we do a joint analysis of the likelihoods from 
dA^z) and those from number counts. To forecast the 
constraints from cluster surveys we use MCMC simula- 
tions in the parameter space of 6 cosmological parameters 
- rim, Wo J Wa, h, Us and fib and 3 scaling parameters - 
A, a and 7. 



IV. RESULTS AND DISCUSSIONS 



Our main results are summarized in Tables |T] and |TT1 
which list the 1-a uncertainties on the cosmological pa- 
rameters. The numbers in bold highlight the improve- 
ment in the dark energy constraints on adding dA{z) to 
^ . The benefits of having additional information from 
dA{z) is also illustrated in figures [5] and [3] In all cases, 
we consider eROSITA as the X-Ray survey and either 
Planck or ACT/SPT as the SZE survey. 




FIG. 2: Cosmological constraints from ACT/SPT: The 1-cr 
regions in the order of overlap (decreasing area) are, number 
counts - with mass follow-up(no priors on A and a) with 
AM ~ 30 - 100% (red) and AM - 15 - 50% (dark-blue); 
only number counts (light-green); with dxiz) (TVbcam = 2, 
578 clusters) (dark-green); with CMB priors on Q.m and as 
(orange); with 557 SNe from SCP Union 2 dataset (light- 
blue). 



A. Constraints from d^iz) added to ACT/SPT ^ 

In figure [21 we see a significant improvement in cosmo- 
logical constraints from ACT/SPT, when dA for A^boam = 

2 is added to ^ . As seen in Table [Tj the 1-a region in 
the dark energy plane WQ-Wa shrinks by factors of 3.8, 
1.7 and 1.6 when the dA{z) datasets with A'bcam = 2, 3 
and 4 are added to ^ ; while in the r2,„ - wq plane it de- 
creases by factors of 4.1, 1.8 and 1.1. The corresponding 
marginalized constraints on a single equation of state w 
improves from 0.20 to 0.15 and 0.18 for A'bcam = 2 and 

3 respectively when dA{z) is added to ^ . There is not 
much improvement with the smallest dataset of Abeam — 
4, since the dA{z) sample has only 25 clusters now. 

Since SNe la are also used to measure distances 
(d-L^z))^ we compare our results with the benefits of 
adding the Union 2 compilation of the Supernova Cos- 
mology Project (SCP) dataset in column 6. As expected, 
with the SNe data of a similar size but with much smaller 
errors (7-10%), we get tighter constraints on dark en- 
ergy with Awo = 0.15 and Awq = 1.0. Column 7 of the 



TABLE I: Comparison of 1 — ct parameter constraints from ACT/SPT ^ and its overlap with eROSlTA. Nd is the number 
of clusters for which dA{z) is measured and NsNe is the number of SNe. Columns 3-5 hsts improvements in the constraints 
obtained from ACT/SPT ^ when dA{z) from various datasets are added. Column 6 shows the constrains when CMB priors 
are imposed on flm and as from the WMAP7 results. The last two columns show the constraints when there are no priors on 
the scaling parameters A and a, but a mass follow-up of 97 clusters is added with (randomly distributed) errors on the masses 
of 15-50% and 30-100% for the follow-up's 1 and 2 respectively. 



Parameter only ^ i§+dA '^+dA ^+dA ^+SNe 


^-hCMB priors ^-hfoUow-l ^-ffollow-2 


Nci/NsNe 995 578 118 25 557 


97 97 


(iVbeam = 2) (A^beam = 3) (iVbeam = 4) 


on Q.m and as AM ~ 15 - 50% AM ~ 30 - 100% 



An„ 


0.067 


0.036 


0.053 


0.052 


0.029 


0.018 


0.073 


0.090 


Awo 


0.79 


0.36 


0.57 


0.63 


0.15 


0.40 


0.81 


0.80 


AWa 


3.5 


2.0 


2.8 


2.7 


1.0 


2.0 


3.4 


3.3 


Acrg 


0.12 


0.10 


0.10 


0.11 


0.10 


0.05 


0.06 


0.07 



TABLE II: Comparison of 1 — cr parameter constraints from 
eROSITA ^ & its overlap with Planck. 



Parameter 

Nd 


oniy^ 

1994 


^+rfA 
l^-i *bcam 

1829 


2) (A^bcam = 

951 


'-t+dA 
3) (Afbcam = 4) 

202 


An,„ 
Alt) 
Ao-g 


0.027 
0.17 

0.086 


0.027 
0.10 

0.099 


0.026 
0.13 

0.083 


0.028 
0.15 

0.093 



same table lists the effect of using CMB constraints from 
the WMAP7 results on ^^(=0.019) and cr8(=0.059) with 
the ^ data from ACT/SPT to further improve the con- 
straints on dark energy with Awg = 0.4 and Awa = 2.0. 
These constraints are very similar to the ones through 
dx{z) for iVboam — 2; compare column 3 and 7 in Table 
III 

Next, we compare our constraints with those that 
would be obtained from dedicated follow-up observations 
of cluster masses requiring detailed X-Ray, SZE or galaxy 
spectroscopic observations [16]. Adding a mass follow-up 
to ^ constraints the cluster scaling relation and is simi- 
lar to putting priors on the scaling parameters. We build 
two mock catalogs of 97 clusters each, one with errors 
between 15-50% and the other with 30-100% errors. We 
find that the constrains from both the mass follow-up's 
(now without priors on the scaling parameters A and a) 
are very similar to those obtained using priors of the scal- 
ing relations, compare column 2 with 8 and 9 in Table 
|T1 This also implies that addition of dA{z) can further 
improve constraints even when a mass follow-up is done. 
Compare columns 3-5 with 8 and 9. 



B. Constraints from dA{z) added to Planck 



dN 
dz 



The eROSITA survey will detect all the clusters seen 
by Planck and hence together, these surveys will yield a 
large number of candidates for estimating dA{z). In our 
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FIG. 3: Constraints obtained from Planck ^ added to dAiz) 
measured from the overlap of Planck plus eROSITA survey 
areas. The figures show 1-a regions in the order of overlap 
(decreasing area) as follows: just ^ (yellow); ^ with c(a(z) 
for A''beam = 4 (light-green), 3 (indigo) and 2 (magenta) re- 
spectively. 



mock catalog the number of cluster datasets for which 
dA can be measured, for this combination, are found to 
contain 1829, 951 and 202 halos for iVbcam = 2, 3 and 
4, respectively. However, in contrast with ACT/SPT, 
Planck would be able to detect only the most massive 
clusters due to a lower resolution and mostly the ones 
occurring at lower redshifts. There would be a only very 
few clusters at higher redshifts (> 0.6) and hence we 
consider the constraints only on a constant equation of 
state w. Table In] lists the results for the improvement in 
the constraints from ^ when dAiz) from an overlap of 
Planck and eROSITA are added. We see here that Aw 
decreases by factors of 1.7, 1.3 and 1.1 with the addition 
oi dA{z) datasets with A^bcam = 2, 3 and 4 respectively. 



V. CONCLUSIONS 

We show that just using number count observa- 
tions from the ongoing SZE galaxy cluster surveys like 
ACT/SPT would not be sufficient to put tight constraints 
on the cosmological parameters, especially the dark en- 
ergy parameters. However these constraints can be sig- 
nificantly improved by doing a joint analysis with dA{z) . 
These dA{z) can be constructed out of clusters detected 
jointly in X-Ray and SZE observations having overlap- 
ping areas. The optical follow-up providing the redshifts 
for ^ will, also, naturally provide the redshifts for dA{z) 
. With the current and upcoming large area surveys one 
will be able to get thousands of clusters providing us with 
us with dA at various redshifts without much effort. 

We find that adding dA{z) to number count observa- 



tions always improves the dark energy constraints, from 
^ alone by factors of 1.5 to 4- This leads to better 
constraints not only on dark energy but also on the pa- 
rameter flm- Even when a targeted mass follow-up of 
clusters helps in breaking the cluster-cosmology degen- 
eracies, addition of dA{z) helps in further tightening of 
the cosmological constraints. Moreover addition of dA{z) 
improves dark energy constraints comparable to the im- 
provement brought by adding CMB priors on as and fi™. 
Thus, our proposal of adding the dA{z) data to cluster 
number counts provides a natural way of improving the 
cosmological constraints using clusters alone. 
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